X-ray Observations of High-B Radio Pulsars 



S. A. Olauseni, W. W. ZW, J. K. Vogel^, V. M. Kaspi^ A. G. Lyne"^, C. M. Espinoza"^, B. W. 
Stappers^, R. N. Manchester^, & M. A. McLaughlin'^ 

ABSTRACT 

The study of high-magnetic-field pulsars is important for examining the relationships 
between radio pulsars, magnetars, and X-ray-isolated neutron stars (XINSs). Here we 
report on X-ray observations of three such high-magnetic-field radio pulsars. We first 
present the results of a deep XMM-Newton observation of PSR J1734— 3333, taken 
to follow up on its initial detection in 2009. The pulsar's spectrum is well fit by a 
blackbody with a temperature of 300ib60 eV, with bolometric luminosity Lbb = 2.0^q7 x 
2q32 gj-gs^^ ~ 0.0036E' for a distance of 6.1 kpc. We detect no X-ray pulsations from the 
j source, setting a la upper limit on the pulsed fraction of 60% in the 0.5-3 keV band. We 

compare PSR J1734— 3333 to other rotation-powered pulsars of similar age and find that 
it is significantly hotter, supporting the hypothesis that the magnetic field affects the 
observed thermal properties of pulsars. We also report on XMM-Newton and Chandra 
observations of PSRs B1845— 19 and JlOOl— 5939. We do not detect either pulsar, 
setting 3a upper limits on their blackbody temperatures of 48 and 56 eV, respectively. 
Despite the similarities in rotational properties, these sources are significantly cooler 
than all but one of the XINSs, which we attribute to the two groups having been born 
with different magnetic fields and hence evolving differently. 

Subject headings: pulsars: general — pulsars: individual (PSR B1845— 19, PSR JlOOl— 5939, 
PSR J1734— 3333) — stars: neutron — X-rays: stars 
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1. Introduction 

Over the past few decades, the boom in X-ray and gamma-ray astronomy has led to a significant 
increase in our knowledge about the neutron star family. Previously, the only kind of isolated 
neutron stars known were radio pulsars. The latter are also known as 'rotation-powered pulsars' 
(RPPs) because their luminosities are generally much lower than E, the rate of energy loss due to 
rotation. Since then, X-ray observations have led to the discovery of several new clas ses of i s olated 



neutron stars, including magnetars and X-ray-isolated neutron stars (XINsfl); see iKasoil jlOld ) 



for a review. All of these classes exhibit properties different from those of conventional RPPs. 

Magnetars, including the soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs), 
are isolated, slowly rotating neutron stars. The known object^ have spin periods in the range 
of 2-12 s, X-ray luminosities that are in many cases much higher than their spin-down lumi- 
nosities, and, assuming conventional magnetic dipole braking, typically very high magnetic fields 
{B ~ 10^^-10^^ G) . It is ger ierally believed that their X - ray luminosities are fueled by the decay of 
these high magnetic fields (iThompson &: DuncanI Il995l . Il996l : [Thompson et al.l 12002) ■ Magnetars 
also show great X-ray variability and activity, ranging from short SGR-like X-ray/7-ray bursts to 
major, slow-rise and long-lived X-ray flares, that are thought to be produced by fracturing of the 
crust and reconfiguration o f t he magnetic field lines. For reviews of magnetars, see for example 
Woods fc Thompson! (|2006l ) or iMereghettil ^20oi ). 



The XINSs (see iHaberll 120071 and iTuroUal l2009l for reviews) are a group of nearby (distance 



< 500 pc) neutron stars characterized by soft thermal X-ray spectra and no detected radio emission. 
Those XINSs with detailed timing measurements have spin periods in the ran ge of 3-11 s, high in- 



ferred magnetic fields (B ^ (1-3) x 10^^ G ) , and characteristic ages of ~10^ yr ( Kaplan &: van Kerkwijk 



2009 , 2011 ; van Kerkwijk Sz Kaplan 20081 ) . T he XINSs are considerably hot ter and more X-ray lu- 
minous than similarly aged ordinary RPPs ([Kaplan &: van Kerkwijkl l2009l ). although they show 
no bursting activity and are less luminous than the magneta.rs. T herefore, it has been proposed 
(|Kaplan van KerkwijM bood : IPons et al.l bood : IPopov et all boid ) that the coohng of XINSs is 
affected by heating of the neutron-star surface due to magnetic- field decay. If correct, this suggests 
that XINSs may be evolutionarily related to the magnetars; one possibility is that XINSs are mag- 
netar descendants, being older versions of neutron stars that were born with 10^^-10^^ G magnetic 
fields which have slowly decayed with time. 

The hint of magnetic field decay heating in XINSs suggests that similar processes may be 
occur ing in radio pulsars that have ages and magn etic fields comparable t o those of the INS s. 



Thanks to radio surveys over the past decade (e.g. [Manchester et al.l I2OOII : iMorris et al.l |2002| ) 



a class of radio pulsar having inferred magnetic fields comparable to XINSs and even bona fide 
magnetars has emerged. What physically differentiates these high-i? RPPs from these other classes 



^Also known as X-ray dim isolated neutron stars (XDINSs). 

^See the McGill SGR/AXP Online Catalog at jhttp: 7/wwwTphysics .mcgill . ca/-pulsar/magiietar /main. html] 
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is unclear; aside from their radio emission, which was selected for via their discovery in radio surveys, 
their X-ray emission properties may well be similar to those of the XINSs and magnetars. Indeed, 
in 2009 the High Time Resolution Universe survey discovered a radio pulsar whose rotational and 
X-ray properties identified it as a magnetar ([Levin et al.ll2010l ). 



Here we present X-ray observations of three high-i? radio pulsars: PSRs J1734— 3333, B1845— 19, 
and JlOOl— 5939. These as a group are key sources for X-ray study, as they share spin, and hence 
spin-inferred properties with the XINSs and magnetars. Thus they might be expected to share 
X-ray properties, too. 

PSR J1734— 3333 has period P = 1.1 7 s and period der ivative P = 2.3 x 10~^^ and was 
discovered in the Parkes Multibeam Survey (jMorris et al.ll2002l ). Its spin parameters imply a spin- 



down luminosity oi E = 4tt IP/P = 5.6 x 10 ergs" , characteristic age of r = P/2P = 8.1 kyr, 
and an inferred surface dipolar magnetic field of = 3.2 x 10^^(i-*P)^^^G = 5.2 x 10^'^ G, which 
is among the highest of all known radio pulsars and similar to those of bona fide magnetars such 
as IE 2259-F586 (B = 5 . 9 x 10^=^ G: jcaspi et al.l [l999 ) and Swift J182 2.3-1606 {B = 4.7 x 10^^ q. 
Livingstone et al.l l201ll : IScholz et al.l |2012| : but see iRea et al.l |2012| ) . It has a radio dispersion 
measure (DM) of 578 p c cm~^ which, based on the NE2001 model for Galactic free electron density 
( Gordes Sz Lazid l2002l ). gives a best-estimated distance to the pulsar of 6.1 kpc (although these 



distance estimates typically h ave large uncertainties of 25% or more). Based on its unusually low 
braking index (n = 0.9 ± 0.2), lEspinoza et al.l (j201ll ) suggested that the pulsar's magnetic field is 



growing, i.e. its trajectory on a conventional P/P diagram is up and to the right, towards the 
region occupied by the magnetars. To search for anomalous X-ray emission from the pulsar, a 
short XMM-Newton observation was taken in 2009. The source was detected, and although there 
were indications that it was anomalously hot, it was too fa int and the observat ion too short to 
determine the spectrum unambiguously or detect pulsations (lOlausen et al.ll2010l ). 



PSR B1845— 19 is a 4.31-s period radio pulsar discovered in the Molonglo pulsar survey 
(|Manchester et al.lll97d ). It has P = 2.33 x 10~^^, ^ = 1.1 x 103iergs~\ r = 2.9 Mvr, and 
B = 1.0 X 10^3 G. It has a DM of 18 pccm , implying a distance of 0.75 kpc (jGordes &: Lazid 



2OO2I ). The pulsar's magnetic field, age, E and distance are close to those of the XINSs. 



PSR JlOOl— 593 9 was another pulsar discovered in the Parkes Multibeam Pulsar Survey 
( Lorimer et al.l 120061 ). with period P = 7.73 s and period derivative P = 5.99 x 10~^^. As 
with PSR B1845-19, the low E (5.1 x 10^°ergs"^), age (r = 2.1 Myr), and high magnetic field 
{B = 2.2 X 10^3 G) of PSR JlOOl-5939 are similar to the properties of the known XINSs, although 
at a distance of 2.7 kpc (as determined from its DM of 113pccm~^) it is farther away than any of 
them. 

In this paper, we report on X-ray observations of these three high-i? RPPs, done with the 
aim of searching for magnetar- or XINS-like emission as might be expected if magnetic-field-decay 
heating is occuring. Section [2] describes our X-ray observations. Section [3] our results, and we 
discuss the implications of our findings in Section HI 
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Observations 



2.1. PSR J1734-3333 



A deep, 125-ks observat ion of PSR J1734- 3333 was carried out on 2011 March 11-12 with the 
XMM-Newton observatory ( Jansen et al.|[200ll ). We also reanalyzed t he EPIC pn data from the 
10-ks 2009 XMM-Newton observation of the source (jOIausen et al.ll2010l ). but not the MOS data, as 
they had neither the proper time resolution nor sufficient sourc e counts to meaning fully contribute 
to our analysis. In both observations, the EPIC p n camera dstriider et all boOlh was operating 



in large-window mode and the EPIC MOS cameras (jTurner et al 



2001 



in full-window mode; the 
medium filter was in use for all three cameras. Details of the two observations, as well as those 
described in Sections 12.21 and \'2.'6\ are summarized in Table [TJ The data from both observations 
were analyzed with the XMM Science Analysis System (SAS) version ii.o.dl. To search for times 
of high background flaring that are known to sometimes affect XMM-Newton data, we extracted 
light curves from over the entire field of view of all three cameras. The 2011 observation was 
heavily affected, with over half the exposure length contaminated by background flares. The 2009 
observation showed no such problems. 



2.2. PSR B1845-19 

PSR B1845-19 was observed for 38 ks by XMM-Newton on 2011 March 16. The pn camera 
was operating in large-window mode with the thin filter, and the MOS cameras were operating in 
full-window mode also using the thin filter (see Table [1] for further details) . The data were reduced 
with XMM SAS 11.0.0, and times of high background flaring were removed. Because the two MOS 
cameras are similar instruments, we combined their images for the analysis of these data. 



2.3. PSR JlOOl-5939 



On 2011 October 8, the Chandra X-ray Observatory observed the position of PSR JlOOl— 5939 
(see Table d]) fo r 19 ks, using the ba ck-illuminated S3 chip of the Advanced CCD Imaging Spec- 
trometer (ACIS: lGarmire et al.ll2003l ) as aimpoint. The 1/8 subarray mode was chosen to achieve a 
time resolution of 0.44 s, sufficient for timing a pulsar of period 7.73 s. For our analysis, we worked 
with the level 2 event files, which were created from level 1 event lists by including good time 
intervals as well as information about issues such as bad pixels and cosmic rays. The data were 
analyzed using the tools provided in CIAO0 4.3. 



See [http : //xmm. esac . esa. int/sas/| 
^See jhttp: //cxc .hajTvard. edu/ciao/| 
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3. Results 
3.1. PSR J1734-3333 

3.1.1. Imaging 

Using the SAS tool edetect_chain to perform a blind search for point sources, we detected 
the X-ray counterpart of PSR J1734— 3333 in both the 2009 and 2011 observations, with the re- 
ported count rates in both observations being consistent with each other. Our best-fit position 
for the X-ray source is the one reported by edetect_chain for the 2011 observation: (J2000) 
R.A. = 17''34™27!03(4), decl. = -33°33'22'/4(5), where the hsted uncertainties are the statisti- 
cal errors only. It is consistent with the radio tir ning position o f R.A. = 17'^34™26!9(2), decl. = 



33°33'20"(10) (jEspinoza et al.ll201ll ). However, lOlausen et al.l (|2010l ) also identified an optical 



source (NOMAD Catalog] ID 0564-0621454) that lay within the 2009 X-ray error circle and de- 
termined that either it or the radio pulsar could be associated with the X-ray source, but not 
both. We find that the updated X-ray position excludes the optical source, even when including^ 
XMM -Newton^s 2" absolute pointing uncertaintjjf], strengthening the conclusion from Olausen et al 



(j20ld ) that the X-ray source is the counterpart to PSR J1734— 3333. We also found that the radial 
profile of the X-ray source is consistent with the XMM point-spread function (PSF). Therefore 
there is no evidence of extended emission. 



3.1.2. Timing Analysis 

To search for X-ray pulsations from PSR J1734— 3333, we extracted counts in the 0.5-3 keV 
energy band from a region of 30" radius centered on the source. The MOS cameras, operating with 
a time resolution of 2.6 s, were unsuitable for this analysis, so only the pn data were used. Within 
the 2011 source region we found 576 counts, and by extracting counts in the same energy range 
from source-free regions on the same CCDs, we determined that 352 it 12 of those counts were from 
the background. Following the same procedure for the 2009 data we determined that of the 150 
total counts in the source region, 100 it 6 were from background photons. After barycentering the 
source region events with the SAS tool barycen, we folded them into eight phas e bins based on 



a rad io timing ephemeris obtained using the Jodrell Bank 76-m Lovell telescope (jEspinoza et al 



2011 



). Since the time span covered by the radio ephemeris included the epochs of both the 2009 



and 2011 observations, we could produce a single summed light curve from both sets of data. 

Fitting the folded curve to a straight line gave a best-fit of 11.1 for 7 degrees of freedom, 
or a 13% chance that the folded curve could result with no signal in the data. Therefore, no 



' ^http : //www . usno . navy . mil/USNO/ astrometry/opt ical- IR-prod/nomad] 
^See |http: //xiimi2 . esac . esa. int /docs/document s/CAL-TM- 0018 .pdf | 
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significant pulsations were detected. We also performed additional searches, attempting to improve 
the signal-to-noise ratio by using an energy range of 1-2 keV or a source region of 20" radius, as 
well as searching around the predicted period, but in all cases there were again no detections of 
pulsations. 

To find an upper limit for the pulsed fraction, we simulated event lists with the same number 
of total counts as found in the source region. The simulated signal had a sinusoidal profile with 
a random phase and had a user-specified area pulsed fraction, where the latter is defined as the 
ratio of the pulsed part of the profile to the entire profile. We found that a signal with area pulsed 
fraction of ~0.23 would be detected with >3a significance 68% of the time. Therefore, because 
62% of the counts in the source region are from the background, we estimate that the la upper 
limit on the area pulsed fraction of PSR J1734— 3333 is 0.6 in the 0.5-3 keV energy range. 



3.1.3. Spectral Analysis 



We extracted the spectrum of PSR J1734— 3333 from the 2011 pn and MOS data using a 
region of 30" radius centered on the X-ray source position. Background regions were chosen from 
nearby, source-free areas on the same CCD as the pulsar, and response and ancillary response files 
were generated with the SAS tasks rmf gen and arf gen. The pn and MOS spectra were grouped to 
have a minimum of 25 counts per bin. We also extracted a spectrum fr o m the 2009 pn data using 
the same source and background regions as described in lOlausen et al.l (j2010l ). grouped to have a 
minimum of 20 counts per bin. 

These four spectra wer e then jointly fit in XSPEC 12.7 .0 using simple power-law and blackbody 
models, as well as with nsa (jPavlov et al.lll995l ) and nsmax (jHo et al.ll2008l ) neutron-star atmosphere 
models, all using phabs to model interstellar absorption. The absorbed power-law model provided 
a statistically acceptable fit to the data (x^ = 34.6 for 44 degrees of freedom), but as it gave an 
unphysically steep photon index of 5.2, we do not consider it any further. Best-fit parameters for 
the other three models are shown in Table [2j 

The best-fit blackbody model is plotted in Figure [H We found a best-fiJll temperature of kT = 
0.30 ± 0.06 keV, corresponding to a column density of A'^h = 0.671q24 x 10^^ cm^^. The blackbody 



radius is i?bb = 0-45!to'2o '^6.1 km, giving a bolometric luminosity of Lbb = 2.0 



cm ■ 

+2.2 
-0.7 



10^2 erg s-\ 

where dg.i is the distance to the pulsar in units of 6.1 kpc. In order to better explore the confidence 
range of kT and Ah for the blackbody model, we plotted their confidence contours in Figure [2j 
The plot shows that the 3a lower limit on kT is 0.18 keV and that lower values of kT require higher 
values of A^h- We note that the maximum column density along the line of sight to the pulsajfl 
is A^H = (1.1-1.4) X 10^2 cm~2, and because the DM distance places the pulsar less than halfway 



^AU errors in this section are 90% confidence intervals. 



Ihttp : //heasarc . nasa . gov/cgi-bin/Tools/w3iih/w3nh ■ pi 
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through the Galaxy, the expected column density is significantly lower than that. As a result, 
models with kT < 0.2 keV which require iVn > 1 x 10^^ cm~^ are further disfavored. 

Both of the neutron-star atmosphere models provide good fits to the data and yield much 
lower best-fit temperatures than the blackbody model, as is typical. The nsa model gives kT°° = 
0.14 lb 0.05 kcV, but because it assumes that the emission is from the entire surface of the star, 
it implies an unreasonably large distance to the pulsar of ~27kpc. Fixing the distance at 6.1 kpc 
results in an even lower temperature of 90 ± 3eV (x^ = 37.6 for 45 degrees of freedom), but 
the required column density rises to Nn = 1.21q"^ x lO^^cm"^, which is not favored as described 
above. Given these difficulties, we also fit the spectrum to the nsmax model for a highly magnetized 
neutron star atmosphere. This model gives fcT^ = 0.13+2;|]^ keV, A^h = 0.88^^;^^ x 1022cm-2, and 
an emission radius of R = 3.2]']2 q '^6.1 km. Finally, we note that the nsa and nsmax models were 
fit using the highest available values for the magnetic field (in XSPEC), 10^^ and 2 x 10^^ G, 
respectively. Since these values are below the 5.2 x 10^^ G field of PSR J1734— 3333, the results are 
not completely reliable. 

3.2. PSR B1845-19 

We employed the SAS tool edetect_chain to search for the X-ray counterpart of PSR B1845— 19, 
but detected no X-ray sources within 1.5' of the pulsar's radio timing position. From the pn image, 
we found 561 counts in a circular region of radius 28" centered on the pulsar position and 781 counts 
from an annular region of outer radius 45" around the source region. From the combined MOS 
image, we found 132 counts in a 27" radius source region and 192 counts in a surrounding annulus 
with outer radius 45". In both cases the source regions were chosen such that they encircle ~80% 
of the energy from the pulsar according to the on-axis PSFs of the instruments. The background 
regions are expected to contain a further --^10% of the source energy. To estimate the count rate 
upper limit of the pulsar, we input these values into the aprates tool from CIAO 4.4 and obtained 
a 3a (99.7% confidence) upper limit of 0.0048 counts s^^ for the pn camera and 0.0029 counts s~^ 
for the combined MOS camera. All counts and count rates are for the 0.2-3 keV energy band. 

In order to convert the count rate upper limit into an upper limit on the pulsar's surface 
temperature kT, we assume the pulsar's spectrum can be described by an absorbed blackbody 
model. We use XSPEC to find the spectrum for different values of the model parameters (blackbody 
radius i^bb, distance, and Nu) and, based on response and ancillary response files generated for the 
source regions described previously, calculate the predicted count rates for the pn and MOS. We 
select 5, 10, and 15 km as our possible values for i?bbi and we use the DM distance to the pulsar of 
0.75 kpc. The total column density along the pulsar's line of sight^ is Ah = (0.14-0.18) x 10^^ cm"^, 
so we evaluate kT upper limits for a range of column densities with Ocm"^ < Ah < 0.2 x 10^^ cm"^ 
and calculate bolometric luminosity based on kT and i?bb- 

We found that the upper limits from the combined MOS cameras provide less strict constraints 
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than those of the pn camera, so we present only the upper hmits from the pn image of PSR 
B1845— 19 here. These upper hmits, in both kT and bolometric luminosity, are shown in Figure [3j 
Taking the rather conservative values of A^h = 0.2 x 10^^ cm~^ and i2bb = 10 km, our 3a upper 
limits on kT and bolometric luminosity are 48 eV and 6.8 x lO^^ ergs'^ « 6.2E. 



3.3. PSR JlOOl-5939 



No X-ray signal was observed at the expected radio position of PSR JlOOl— 5939. We extracted 
counts over the 0.3-3 keV energy range from a 5" radius circular region centered on the pulsar 
position and 141 counts in the same energy range from a surrounding annular region with outer 
radius 60". The radius of the source region was chosen such that it contains ~99% of the expected 
signal based on the Chandra on-axis PSF; the rest of the potential signal should be present in 
the background annulus. From these values the CIAO aprates tool gives a 3a upper limit on the 
pulsar's 0.3-3 keV count rate of 0.34 x 10~^ counts s~^. 

Using the same procedure as outlined above for PSR B1845— 19, we can obtain upper limits 
for the surface temperature of PSR JlOOl— 5939 from our count rate upper limits. We again choose 
blackbody radii of -Rbb = 5, 10, and 15 km, and use the DM distance to the pulsar of 2.7 kpc. 
Based on the DM distance, we also estimate the column density to be i Vn = 0.33 x 10^^ cm^^ 



from the 3-dimensional Gala ctic extinction (Ay) mode l of iDrimmel et al.l (j2003l ) and the relation 
Nn/A]/ = 1.79 X 10^^ cm~^ in lPredehl Sz SchmittI ( 19951 ). Therefore, we evaluate kT and luminosity 
upper limits for 0.1 x 10^^ cm~^ < iVn < 0.55 x 10^^ cm~^. Note that the total A^'h along the hue 
of sight to the pulsar^ is (0.80-0.85) x 10^^ cm 



,-2 



Figuredlshows our results for the kT and bolometric luminosity upper limits of PSR JlOOl— 5939. 
In particular, for our estimate of A^h = 0.33 x 10^^ cm~'^ and assuming a blackbody radius of 10 km, 
we find 3a upper limits of 56 eV for kT and 1.26 x 10^^ 

25E for the bolometric luminosity. 



4. Discussion 

4.1. PSR J1734-3333 

We have reported on a new long XMM-Newton observation of the young, high-S pulsar PSR 
J1734— 3333. We found a more precise position for the X-ray counterpart that falls within the error 
ellipse of the pulsar's radio timing position. We were unable to detect X-ray pulsations, and the 
la upper limit of 60% (0.5-3 keV) that we set on the pulsed fraction is not very constraining. 

PSR J1734— 3333 has a thermal spectrum, well fit by a blackbody model with an unusually 
high temperature, kT = 0.30 ± 0.06 keV, and corresponding blackbody radius of 0.451q'2o de.i km. 
The bolometric luminosity is Lbb = 2.0 x lO^^dg ^ ergs~^, giving an X-ray efficiency of L^^h/E = 
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0.0036d|i. The pulsar's small emission radius is suggestive of thermal emission from heated polar 
caps, but such models predict X - ray effi ciencies less than 10~^. For example, the polar cap reheating 
model of iHarding Muslimovl (|200lh predicts Lyj^jE « 3 x lO""^ for PSR J1734-3333, an order 
of magnitude below what we observe. Of course, if the DM distance is incorrect and the true 
distance to the pulsar is within 2kpc (de.i < .3), then the X-ray efficiency drops to roughly what 
is predicted for polar cap reheating. However, iHarding &: Muslimovl (l200ll ) warn their model may 
not apply to pulsars with i? > 4 x 10^^ G so these predictions are unreliable. 

The pulsar's blackbody temperature of .30 keV is almos t three times as high as expected 
from a minimal cooling model (0.07-0.11 keV; iPage et al.l 120061 ) given its age. Even the neutron- 
star atmosphere models give best-fit temperatures that are too high, kT^ = 0.13-0.14 keV, though 
their uncertainties do overlap with the range of predicted cooling temperatures. For instance, the 
nsa model with distance fixed at 6.1 kpc gives kT°^ = 0.09 keV. However, note that as with the 
simple blackbody, lower temperatures are correlated with higher iVn in the atmosphere models too, 
which can disfavor that region of parameter space as explained in Section 13.1.31 

Although our ignorance regarding the correct atmosphere model for PSR J1734— 3333 pre- 
cludes meaningful comparisons of its X-ray emission with models of neutron star cooling, we can 
nevertheless compare our results with those of other RPPs of comparable age. We do this in 
Figure [5] where we plot blackbody temperature versus characteristic age for RPPs, XINSs, and 
magnetars. High-i? sources are shown in red and yellow (where "high" means B > 10^'^ G). Al- 
though a blackbody model neglects atmospheric effects and hence is almost certainly incorrect for 
these sources, it provides a simple spectral parameterization by which to describe all source spectra 
and allows a consistent comparison of spectral properties, even if ultimately not fully physical. 
Figure [5] shows that PSR J1734— 3333 has a much higher blackbody temperature than nearly all 
other radio pulsars for which this q uantity has been mea sured, with the possible exception of PSR 
J1119-6127, another high-5 RPP (jGonzalez et allboosl : Ng et al. 2012, in preparation). It is even 
hotter than the magnetar XTE J1810— 197, although its temperature is well below those of the 
other magnetars. In particular, however, it is far hott er than multiple other RPPs of comparable 



age but smaller magnetic field. Indeed, as noted by IZhu et al.l (|201ll ). high-i? RPPs in general 



show a trend toward higher blackbody temperatures when compared with lower-B RPPs of similar 
age. This suggests that the magnetic field affects the observed thermal properties of RPPs, as is 
expected if the star is actively heated by field decay. However, enhanced thermal emission could 



also be a result of pass ive effects such as magnetically altered thermal conductivity (jGeppert et al 



20061 : iPage et al.l 120071 ) without any active field decay. Distinguishing between active and passive 
explanations for the enhanced thermal emission in high-i? RPPs could be achieved via detailed 
modelling, or perhaps by observing other signatures of active field decay, s uch as magnetar-lik e 
bursts in high-S RPPs. Indeed this has been observed in one source already ([Gavriil et al.ll2008l ). 
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4.2. PSRs B1845-19 and JlOOl-5939 



We have reported on X-ray observations of two rotation-powered pulsars, PSRs B1845— 19 and 
JlOOl— 5939, taken with the XMM-Newton and Chandra telescopes, respectively. In both cases we 
failed to detect the pulsars and set upper limits on their kT, assuming a blackbody spectrum. These 
two sources were selected for observation because their spin-derived properties are very similar to 
those of the XINSs, but as shown in Figure EJ this is not the case in the X-ray band. The kT upper 
limits we find for our two pulsars imply that they are significantly cooler than the bulk of the XINS 
population, although we cannot exclude the po ssibility that their t emperatures are similar to that 
of the coolest known XINS, RX J0420.0-5022 (|Haberl et al.l 120041 ). 



In the magneto-thermal evolution model of iPons et al.l (j2009l ) and iPopov et al.l (j20ld ) , the 
XINSs are born with magnetic fields much higher than what we currently measure (^10^"^ G, possi- 
bly including a significant toroidal component). As they age, the field decays and heats the surface 
of the star. A reasonable explanation, then, for PSRs B1845— 19 and JlOOl— 5939 to have much 
lower temperatures than most of the XINSs, is that they have similar ages but were born with 
lower magnetic fields. They would therefore experience less heating from magnetic-field decay and 
their thermal evolution would more closely resemble t he standard cooling curves of ii on-magnetic 
neutron stars. This explanation was also suggested by iKaplan Sz van KerkwiikI (j201ll ) to account 
for the low temperature of RX J0420. 0—5022, and indeed the upper limits on our sources do allow 
for the possibility that they and RX J0420. 0—5022 followed similar evolutionary tracks. 

Another possible explanation is that PSRs B1845— 19 and JlOOl— 5939 do not have similar ages 
to the XINSs. The above models assume t hat the XINSs have true ages of ~0.5-lMyr, based on 
kinematic age estimates for some of them (|Tetzlaff et ahlboilh . These kinematic ages are smaller 
than those inferred from spin-down (1.5-4 Myr). If the kinematic ages are correct, then it is possible 
that our radio pulsar sources have true ages considerably closer to their characteristic ages (possibly 
as a result of different evolution as above). In that case, they have lower temperatures than most or 
all of the XINSs sim ply because they are older and have cooled further. For a related scenario see 
Turolla et al.l (1201111. who, b ased on magneto-thermal evolution, model the low-field magnetar SGR 
0418+5729 (iRea et al.ll2010l ) as being ~l-1.5Myr old and having undergone substantial magnetic- 
field decay. 



5. Conclusion 

In summary, we have presented results from XMM-Newton and Chandra observations of three 
high-S radio pulsars: PSRs J1734-3333, B1845-19, and JlOOl-5939. Our observation of PSR 
J1734— 3333 was taken to follow up on its detection in 2009. Although we do not detect X- 
ray pulsations from the source, we find that it has a thermal spectrum with kT = 300 it 60 eV, 
significantly higher than predicted by cooling models. Atmosphere models (nsa and nsmax) yield 
lower temperatures that are consistent with standard cooling, but the nsa model in particular 
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predicts too high of a distance or column density. Comparing PSR J1734— 3333 with other RPPs, 
we find it has a significantly higher blackbody temperature than nearly every other radio pulsar 
for which such a measurement has been made, particularly other RPPs of similar age but lower 
magnetic field. This result suggests that a high magnetic field can affect the observed thermal 
properties of pulsars, for example by heating of the surface due to magnetic field decay, or by 
passive atmospheric effects. PSR B1845— 19 and PSR JlOOl— 5939 were not detected in our X-ray 
observations. Assuming blackbody emission from a 10-km radius, we derive 3a upper limits on their 
temperatures of 48 and 56 eV, consistent with standard cooling curves. The two radio pulsars have 
similar spin-derived properties to the XINSs, but the temperature limits we find are significantly 
below the observed temperatures of all the XINSs except RX J0420. 0—5022. We conclude that 
these radio pulsars may have been born with lower magnetic fields than any of the XINSs (with 
the possible exception of RX J0420. 0—5022), experienced less magnetic-field-decay heating, and 
therefore their evolution more closely resembled that of ordinary radio pulsars. 

Part of this work was performed under the auspices of the U.S. Department of Energy by 
Lawrence Livermore National Laboratory under Contract No. DE-AC52-07NA27344. VMK holds 
the Lorne Trottier Chair in Astrophysics and Cosmology, and a Canada Research Chair, a Killam 
Research Fellowship, and acknowledges additional support from an NSERC Discovery Grant, from 
FQRNT via le Centre de Recherche Astrophysique du Quebec and the Canadian Institute for 
Advanced Research. Support for this work was provided by NASA through Chandra award GOl- 
12083X. 
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Table 1. Summary of X-ray Observations 



Target 


Telescope 


Obs ID 


Date 


Detector 


Time Resolution (s) 


Exposure (ks)*^ 


J1734-3333 


XMM 


0553850101 


2009 Mar 9 


pn 


0.048 


8.7 




XMM 


0653320101 


2011 Mar 11 


pn 


0.048 


42.8 










MOSl 


2.6 


59.8 










MOSl 


2.6 


64.5 


B1845-19 


XMM 


0653300101 


2011 Mar 16 


pn 


0.048 


28.2 










MOSl 


2.6 


27.3 










M0S2 


2.6 


27.3 


JlOOl-5939 


Chandra 


12561 


2011 Oct 8 


ACIS-S 


0.44 


17.2 



'^The exposure time is dead-time corrected and has intervals of high background flaring removed. 



Table 2. Spectral Models for PSR J1734-3333 



Parameter Blackbody nsa^ nsmax 



Nh (10^2 cm-2) 


fi7+0-35 
^•"'-0.25 




U.OO_Q 29 


0.88lg:« 


kT"^ (keV) 


0.30 ± 0.06 


0.14 ±0.05 


0.131^1 


(km) 


^•^"^-0.20 


13 (fixed) 


q 9+8.9 
o.z_2 Q 


Distance'^ (kpc) 




271^0 




(dof) 


34.7(44) 


34.6(44) 


34.3(44) 
9.8+i-? 


/abs (10-15 erg s-icm-2)d 


9.9 ± 1.2 


9.9lii 


/unabs (10-14 erg s-icm-2) 


d 3.6+n 


5.611? 


6.61^1/ 



^The nsa model for a pulsar with B = 10 G. 

•^The nsmax model for a pulsar with B = 2 x 10^^ G. 

'^The distance to the pulsar is fit for in the nsa model, while in the 
other two models the DM distance of 6.1 kpc is used to determine 
the radius of emission. 



'I Absorbed and unabsorbed fluxes are given in the 0.5-2.0 keV 
band. 
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2011 pn 
2011 M0S1 
2011 M0S2 
2009 pn 



4x10-3- 



2x10" 




-2x1 0"'' - 



Energy (keV) 



Fig. 1.— XMM spectrum of PSR J1734-3333 with the best-fit blackbody model. The bottom 
panel shows the residuals in units of counts s~^ keV~^. Error bars are at the la confidence level. 
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Fig. 2. — la, 2(T, and 3a confidence contours for kT and A^h from fitting an absorbed black- 
body model to the X-ray spectrum of PSR J1734— 3333. The cross marks the best-fit value, and 
the dashed lines mark the total A'h along the line-of-sight to the pulsar according to the Lei- 
den/Argentine/Bonn (labeled LAB) and Dickey & Lockman (labeled DL) Galactic H I surveys. 
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Fig. 3. — 3a upper limits on the blackbody temperature kT (left) and bolometric luminosity (right) 
for PSR B1845-19. 
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Fig. 4. — 3(7 upper limits on the blackbody temperature kT (left) and bolometric luminosity (right) 
for PSR JlOOl-5939. 
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Fig. 5. — Blackbody temperatures vs. characteristic ages for normal pulsars (blue triangles), 
high-S pulsars {B > 10^^ G; red diamonds), XINSs (red squares), and magnetars (yellow stars). 
Magnetars were plotted only if they had a measurement of kT in quiescence with uncertain- 
ties of 0.1 keV or less (IE 1048.1-5937, IE 1547.0-5408, IE 1841-045, IE 2259+586 IRXS 
J170849.0-400910, 4U 0142+61, CXOU J010043.1-721134, CXOU J164710.2-455216, CXOU 
J171405.7-381031, PSR 1622-4950, SGR 1900+14, and XTE J1810-197), based on data taken 
from ij http : //www . physics . mcgill . ca/ "pulsar /magnetar /main . htmlj The data for the other 
sourc es are taken from IZhu et al.l (j201ll ) with the a ddition of PSR J0726— 2612 (jSpeagle et al 



2011 



), an updated tem perature for PSR J1119— 6127 (jNg et al.ll2012l ). an updated timing solution 



for RX J0420.0-5022 (iKaplan fc van Kerkwiik 



2011 



and data from this work. The three high-i? 



pulsars described in this paper are represented by large red diamonds with holes. 



